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Abstract 

Many dynamic systems containing a large number of modes can benefit from adap- 
tive control techniques, which are well suited to applications that have unknown pa- 
rameters and poorly known operating conditions. In this paper, we focus on a direct 
adaptive control approach that has been extended to handle adaptive rejection of per- 
sistent disturbances. We extend this adaptive control theory to accommodate proble- 
matic modal subsystems of a plant that inhibit the adaptive controller by causing the 
open-loop plant to be non-minimum phase. We will modify the adaptive controller 
with a Residual Mode Filter (RMF) to compensate for problematic modal subsystems, 
thereby allowing the system to satisfy the requirements for the adaptive controller to 
have guaranteed convergence and bounded gains. This paper will be divided into two 
parts. Here in Part I we will review the basic adaptive control approach and introduce 
the primary ideas. In Part II, we will present the RMF methodology and complete the 
proofs of all our results. Also, we will apply the above theoretical results to a simple 
flexible structure example to illustrate the behavior with and without the residual 
mode filter. 


INTRODUCTION 


Applications of control theory to flexible aerospace structures have been many 
and varied. The survey [13] provides a foundation for structure control with many 
control approaches and examples. This was based upon a distributed parameter 
approach to control of flexible structures and other very large-scale systems [14]. 
Later work created the idea of a Residual Mode Filter (RMF) to offset the destabi- 
lizing effect of unmodeled modes in a feedback control environment [15] -[17]. 
This RMF-based structure control theory has been applied to the complex control 
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issues for large horizontal-axis utility-sized wind turbines [18]-[21], and is begin- 
ning to be applied to aeronautic problems that currently use notch filters, eg for 
flutter, also we are applying the theory to aircraft control where there are flexi- 
ble modes in the pilot bandwidth, e.g. large civil tilt rotor. 

In this paper, we extend our adaptive control theory [l]-[4], [7] to accommo- 
date modal subsystems of a plant that inhibit the adaptive controller, in particular 
those residual modes that interfere with the almost strict positive real condition. 
The systems we consider will be large dimensioned, linear time invariant ones 
which can be diagonalized or placed into modal form. This will include linear 
flexible structures of many types. Our adaptive Control approach allows for large 
dimensioned systems through a foundational use of Ideal Trajectories so that the 
adaptive controller is of much lower dimension than the plant. 

The modification will use the idea of Residual Mode Filters (RMF) introduced 
for fixed gain controllers in [6]. In this paper the RMF will be used to eliminate 
the effect of modes that prevent the almost strict positive realness of the overall 
system by being non-minimum phase. This is a new use of the RMF idea; in pre- 
vious non-adaptive work the purpose of the RMF was to eliminate or mitigate the 
destabilizing effect of modes unmodeled in the control system design, whereas 
here the RMF is applied to reinstate the minimum phase nature of the plant under 
adaptive control. 

Here in Part I we will review the basic adaptive control approach and introduce 
the primary ideas. In Part II, we will present the RMF methodology and complete 
the proofs of all our results using results from [8]. Also, we will apply the above 
theoretical results to a simple flexible structure example to illustrate the behavior 
with and without the residual mode filter. 


Rejection of Persistent Disturbances 

The Plant used in this theory section of the paper will be modeled by the linear, 
time- invariant, finite-dimensional system: 

jx p =Ax p +B Up +rn D 

|y„ = Cx P ; x p (0) = x 0 

where the plant state X p (f) , is an N p -dimensional vector, the control input vec- 
tor, U p (t), is M-dimensional, and the sensor output vector , y p (t) , is P- 

dimensional. The disturbance input vector, U 0 (l), is Mn-dimensional and will 
be thought to come from the Disturbance Generator : 
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ju d=® Z D 
K =Fz fl ;z fl (0) = z 0 


2 ) 


where the disturbance state, z D {t) , is No-dimensional. All matrices in (l)-(2) 

will have the appropriate compatible dimensions. Such descriptions of persistent 
disturbances were first used in [5] to describe signals of known form but unknown 
amplitude. Equation (2) can be rewritten as in [3] in a form that is not a dynamical 
system, which is sometimes easier to use: 


fu D =0z„ 
W = L &D 


(3) 


where (j) D is a vector composed of the known basis functions for the solution of 
U D =0Z D ,i.e„ </> D are the basis functions which make up the known form of 
the disturbance, and L is a matrix of dimension No x dim( (f) D ). For the analysis 
performed in this paper, the amplitude of the disturbance does not need to be 
known, so (L,0) can be unknown. For a better understanding of the disturbance 
generator, consider the example of a disturbance generator for a step disturbance; 
in the fonn of equation (2), a step disturbance would have 0 = 1 and F = 0 , in 
the form of equation (3), a step disturbance would have <p D = 1 . 

In [5]-[6], as with much ofthe control literature, it is assumed that the plant and 
disturbance generator parameter matrices, (A,B,C,r, 0, F) , are known. This 
knowledge of the plant and its disturbance generator allows the Separation Prin- 
ciple of Linear Control Theory to be invoked to arrive at a State -Estimator based, 
linear controller which can suppress the persistent disturbances via feedback. In 
this paper, we will not assume that the plant and disturbance generator parameter 
matrices, (A,B,C,r,0) , are known. But, we will assume that the disturbance 

generator parameter from (2), F, is known, i.e., the form of the disturbance func- 
tions is known. In many cases, knowledge of F is not a severe restriction, since the 
disturbance function is often of known form but unknown a mplitude. 

Our control objective will be to cause the output of the plant, y (f) , to 

asymptotically track the output of a known reference model, y m (t) . The Refer- 
ence Model is given by 

j K = K*m + X m(°) = X 0 

1 v = C x 

\ y m m m 


( 4 ) 
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where the reference model state, is an N m -dimensional vector. The refer- 

ence model output, y m (t) , must have the same dimension as the plant output, 
y p {t) . The excitation of the reference model is accomplished via the vector, 
u „,(0 , which is generated by 


U = F U : U (0) = Un (5) 

m m m’ / 0 V ' 

It is assumed that the reference model is stable and the model parameters, 
(A„,B .C F„ ) are known. 

As in [5] -[6], we define the Ideal Trajectories for the plant given by (1) as li- 
near combinations of the plant states, the control inputs, and the disturbance in- 
puts: 


Jx» — S n x m + S 12 u m + S , 3 z D 
[u* = S 21 x m + S 22 u m + 


( 6 ) 


where X*(^) is the ideal trajectory, U*(7) is the ideal control, U*(^) and 

fi. = Ax, +Bu t +ru D ; x < ,(0) = x 0 
\ n (7) 

[y* =C x * = y m 

Note that the ideal output, y *(t) , matches the reference model output, 
y m ( t ) . If such ideal trajectories exist, they will produce exact output tracking. 

By substituting the ideal trajectories given in (6) into (7) and by using the dis- 
turbance generator given by (2), the ideal trajectories can be made to match the 
reference model (4) -(5) with the following Model Matching Conditions: 


as; 

+ 

bs; 

= s;a,„ 

as; 

+ 

bs; 2 

=s;B m +s;f,„ 

as; 

+ 

bs; 

+ r© = s;F 

cs; 

= 

c. 


cs; 

= 

0 


cs; 

= 

0 



( 8 ) 
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The model matching conditions given in (8) are necessary and sufficient condi- 
tions for the existence of ideal trajectories. Solutions to these matching conditions 
must exist for later analysis, but explicit solutions need never be known for the 
adaptive controller design. Necessary and sufficient conditions for the existence 
and uniqueness of solutions to (8) are given in [9], We repeat this result here for 
completeness and the proof is given in the Appendix found in Part II. 

Lemma 1: If CB is nonsingular, there exist unique solutions to the Linear 

Matching Conditions (8) when T(s) = C(sI~A) l B shares no transmission 

zeros with the eigenvalues of A m , F , or F . 

The desired control objective is for the output of the plant to asymptotically 
track the output of the reference model. W e define the output error vector as : 

e v = y p ~y (9) 

To achieve the desired control objective, we want e v — — >0. We define 
the state tracking error as follows: 

e, =x p -x» (10) 

Using (7) and (10), we can write the output error vector as: 

e v = y p - y m = y , - y* = Cx p - cx* = Ce* a i) 

Furthermore, if we let All = U — U* , from (1) and (7) we have 

e* = Ae* +BAu (12) 

For analysis purposes, we define a Fixed Gain Controller 

Up = u» + G*e v (13) 

If we use the fixed gain control law (13) in the plant given by (1), combined with 
the definition of X* from (7) and the output error vector in the form of equation 
(1 1), we obtain: 


e,=(A + BG:c)e* 


( 14 ) 


We can summarize the above by the following: 
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Theorem 1: If (A,B,C) is output feedback stabilizable with a gain G e , 

i.e., the eigenvalues of A c = A + bg:c are all to the left of the jco-axis, then 
the fixed gain controller, (13), will produce asymptotic output tracking, i.e., 

e v — - — — >0. 

y t — >oo 

If all the plant parameters, (A,B,C,T,0,F) , are known, then the fixed gain 
controller given by (13) with a state estimator for z D would be adequate for 
asymptotic tracking. Note that output feedback stabilization of (A, B,C) can be 
accomplished when 


M +P + N d > N p (15) 

and (A, B,C) is controllable and observable [9], In (13), detailed knowledge of 
the parameter matrices is not required, suggesting that an adaptive control scheme 
might be possible under our original assumptions that (A,B,C,T,0) are un- 
known and F from (2) is known. 

Consider the plant given by (1) with the disturbance generator given by (3). 
Our control objective for this system will be accomplished by an Adaptive Control 
Law of the form: 


U p - G m X m +G „ U ,« +G e e y +G D ( / ) D ( 16 ) 

where G m ,G u ,G e , and G d are matrices of the appropriate compatible dimen- 
sions, whose definitions will be given later. We develop the gain adaptation laws 
to make asymptotic output tracking possible by first forming the following which 
are intended to simplify our notation: 

ag„=g„-s; 2 

AG m =G -S 21 

1 (17) 

AG e =G e -G* e 

^G d = G d -S 23 L 

The starred gains in (17) are for analysis and come from the ideal trajectory, 
X* , of equation (6) with Z D in the form given in (3), which is then substituted in- 
to the fixed gain controller (13). Using (6), (7), and the adaptive control law (16), 
we can define: 
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Au = u p - u t 

= AG„u m + AG m x m + (AG e + G*)e v + AG D ^ W 
Then, via (11), (12), and (18), with appropriate definitions, we have: 
e t = Ae„ +BAu 

= (a + BG*c)e* + B[AG„ AG m AG e AGJt? ( 19 ) 

= A c e, + BAG 77 

where Tj = [«f„ x l < fof is the vector of available information . We com- 
bine (12) and (T9) to obtain the Tracking Error System: 

fe« = A c e* +BAG/; 

„ ( 2 °) 

[e v = Ce« 

Now we specify the Adaptive Gain Laws : 

G = -e/7 r H (21) 


where H = [h (7 ], i=l,2,...,4 is an arbitrary, positive definite matrix ( i.e. , H> 0). 
This yields 


G„ =-e v u> n 

G m =-e v x > 22 

i • (22) 

G e = e v e > 33 

G D — — c v ^ Z5 h 44 

Our Adaptive Controller is specified by (16) with the above adaptive gain laws 
(22). Note that none of the starred gains used in the earlier analysis appear in the 
realizable control law, (16) and (22). Next we will analyze the stability ofthis con- 
troller. 

The closed-loop adaptive system consists of (l)-(5), (9), (16), and (22). Using 
(20) and (21), we have 



( 23 ) 


e, = A c e* +BAG/; 
AG = G = -e^ r H 
£ y Ce* 


where A c = A + bg‘c We are able to obtain (23) from (21) because 

AG = G — G* where G* =[S 22 S 21 G e S 23 L] is constant (although gen- 
erally unknown). The stability of the nonlinear system (23) can be analyzed using 
Lyapunov Theory . We form the pos itive definite functions : 


1 T 

Fj(e*) = -e*Pe« 

V 2 (AG) = ^trace( AGH AG r ) 


(24) 


where P > 0 is the solution of the following pair of equations: 

Ja^P + PA c = — Q ; Q > 0 

{pB = C r 


These equations are usually known as the Kalman -Yacubovic Conditions . The 
existence of a symmetric positive definite solution of (25) is known to be equiva- 
lent to the following condition: 

T c (,v) = C(.S’I - A r ) ' B is strictly positivereal ( SPR ) (26) 

For a proof of this equivalence, see [12] App. B. The strict positive realness of 
T C (J) means that for some (7 > 0 and for all co real, 

Re T c (-cr + jco) > 0 (27) 

If we calculate the derivatives, V t , along the trajectories of (23), we have, us- 
ing (25), that 
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V x =-^-e[Qe« +e*PBAG// 


1 


= — e*Qe* +e„v 
2 y 

v = AG/; 


(28) 


and 


V 2 =trace( AGH AG r ) 

= trace(-e v ?i T H)(H ! AG r ) (29) 

= -trace(e v v T ) = - trace(e T v v ) = -e [ v 


1 T 

We can form F = V l + V 2 => V = — — e* Qe* with V < 0 . Consequently, 

Lyapunov theory guarantees stability of the zero equilibrium point of (23) and all 
trajectories of (24) will remain bounded. This guarantees that both e* and AG 
are bounded. 

We can summarize the above by the following Closed-Loop Stability Result: 


Tlieorem2: Suppose the following are true: 

(1) AH U m (0 are bounded (i.e., all eigenvalues of F„, are in the closed left- 
half plane); 

(2) The reference model (4) is stable (i.e., all eigenvalues of A„, are in the 
open left-half plane); 

( 3 ) K is bounded ( i.e., all eigenvalues of F are in the closed left-half plane 
and any eigenvalues on the imaginary axis are simple); 

(4) (A, B, C) is Almost Strict Positive Real (ASPR), i.e., 

T c (s) = C(aI — A c ) 1 B is strictly positive real. 

Then e* and AG are bounded and e* — ■ > 0 and 

e = v —v =Ce* >0. 

y J p J m * oo 

See the Appendix in Part II for a proof of Theorem 2. 

This stability analysis shows that asymptotic tracking occurs and the adaptive 
gains remain bounded. It does not prove that AG — — >0. In fact, the gain 

adaptation laws (22) may not converge to the starred gains in (8); however, this is 
not required for the adaptive controller to achieve its goals. 



10 


Conclusions for Part I 


We have reviewed our adaptive control theory here. This theory accounts for 
adaptive model tracking and for leakage of the disturbance term into the Q modes. 
However, the results require that the error system be minimum phase. In Part II, 
we will show how to modify the adaptive control with residual mode filters to deal 
with non-minimum phase systems. 
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Abstract 

In Part n, we extend our adaptive control theory to accommodate problematic 
modal subsystems of a plant that inhibit the adaptive controller by causing the open- 
loop plant to be non-minimum phase. We will modify the adaptive controller with a 
Residual Mode Filter (RMF) to compensate for problematic modal subsystems, there- 
by allowing the system to satisfy the requirements for the adaptive controller to have 
guaranteed convergence and bounded gains. Also, we will apply the above theoretical 
results to a simple flexible structure example to illustrate the behavior with and with- 
out the residual mode filter. 


INTRODUCTION 


In Part II, we continue the development of the adaptive control approach. We will 
keep the consecutive equation numbering from Part I as well as the same reference 
list. We modify the adaptive control using the idea of Residual Mode Filters 
(RMF) introduced for fixed gain controllers in [6]. In this paper the RMF will be 
used to eliminate the effect of modes that prevent the almost strict positive real- 
ness (ASPR) of the overall system by being non -minimum phase. We have main- 
tained the same reference list here as in the previous paper to make both papere 
more readable. 
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Residual Mode Filter Augmentation of Adaptive Controller 


In some cases the plant in (1) does not satisfy the requirements of ASPR. Instead, 
there maybe be a modal subsystem that inhibits this property. This section will 
present new results for our adaptive control theory. We will modify the adaptive 
controller with a Residual Mode Filter (RMF) to compensate for the troublesome 
modal subsystem, or the Q modes, as was done in [6] for fixed gain non-adaptive 
controllers. Here we present the theory for adaptive controllers modified by 
RMFs. In a previous paper, we examined the RMF with adaptive control, but as- 
sumed that there was no leakage of the disturbance into the Q modes [7]. Here we 
will deal with the issue of disturbances propagating through these modes. 

Let us assume that (1) can be partitioned into the following modal form: 


X 


~ A 

0 " 

X 

+ 

“ B~ 

u + 

“ r “ 

_Xq_ 


0 

to 

1 

x Q _ 


1 

O) 

1 

p 

i 

Co 

i 


r i 

fxl 

[C c e\ 

1 

X 

to 

1 


(28) 


Define X p 


X 

, A„ = 

A 

0 

, B = 

B 

, r = 

r 

, C T = 

~ C~ 


5 P 

0 

A Q- 

’ p 

- B Q- 

’ p 

_ sT q_ 

’ p 

c ’ 
L c eJ 


Disturbance Generator 


I Z D — Fz D 

K = dz D 


or z D —L(j) D as before in (2)-(3). 


The Output Tracking Error and control objective remain as in (4)-(5), i.e. 


e v =y p -y„ 


aO. 


However, now we will only assume that the subsystem ( A,B,C ) is Almost 
Strictly Positive Real, rather than the full un -partitioned plant (A p ,B p , C p ), and 
the modal subsystem ( Aq,B q ,Cq ) will be known and open-loop stable, i.e., A 0 


is stable. Also note that this subsystem is directly affected by the disturbance in- 
put. Recall that ASPR means CB> 0 and P(s) — C(sl — A) X B is minimum 
phase. So, in summary, the actual plant has an ASPR subsystem and a known 
modal subsystem that is stable but inhibits the property of ASPR for the full plant. 
Hence, this modal subsystem must be compensated or filtered away. 

W e define the Residual Mode Filter (RMF): 


ixg - A q X q + B Q U p 
be = CqXq 


( 29 ) 
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And the compensated tracking error : 

S y= e y~yQ 

Now we let e Q = Xg> — x^ and obtain: 


^o e Q ^ Q U D 


( 30 ) 


(31) 


Consequently, 

e y = e y -y Q = CAx + C Q x Q - [C Q x Q + C Q e Q ] 
= CAx - CqCq 


(32) 


As in [l]-[2], we define the Ideal Trajectories, but only for the ASPR Subsys- 
tem: 


fx* = Ax* + Su* + r U D 
{ y* = Cx* = 0 


(33) 


with 


{ x* — iS) z q 
z/* — 


This is equivalent to the Matching Conditions: 


\s*F=AS*+Bsl + re 

Icsf = 0 


(34) 


which are known to be uniquely solvable when CB is nonsingular. However, we 
do not need to know the actual solutions for this adaptive control approach. 


Ax = x-x* 


Let 


(A u = u p -u. 


. Then we have 


[ Av = e y = CAx - C Q e Q 


( Ax = A Ax + BAu 
1 Ay = CAx - C Q e Q 


( 35 ) 


because, from (33), y * = 0. This system can be rewritten: 
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A Q_ 
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Au = A 
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_ e Q_ 


+ B \l I + \ .ll ^ 
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Ax 
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-CqI 

e Q_ 
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_ e Q_ 


(36) 


Now we have the following: 

Lemma: 


f 

A 

0 " 


B 


A = 

V 

0 

A o_ 

,B = 

0_ 

w 

1 

U_, 

ll 

u 


is ASPR if and only if 


(A,B,C) is ASPR. 

Proof: CB-^C —Cq ] 


= CB > 0 and 


P(s) = C(sI -A) l B 

(si -Ay 


=[c - c e] 


0 


(sI-A Q y\ 


= C(sl - A)~ l B = P(s) 

is minimumphase. Ehdofproof. 

So there exists G e such that (A c =A +B G e C ,B ,C) is Strictly Positive Real 
(SPR) when ( A,B,C ) is ASPR. Consequently, as is well known from the Kal- 
man-Yacubovic Theorem, there exists P,Q > 0 such that 


| A'P+PA C =-Q 

1 pb = c t 


(37) 


We now write the modified adaptive control law with RMF: 

Up = G u u m + G m x m + G e e y + G D (j) D 

e y =y P -y Q 

< . 

x Q = A q x q + B Q u p 

y<2 - l 'o x q 


( 38 ) 


with modified adaptive gains given by 
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G u = 
G m = 
G e = 

G n = 


—e u h \h >0 

y m u) u 

-e y x m h m ;h m > 0 
-eXh e ;h e > 0 


( 39 ) 


Finally, we have the following stability result: 

Theorem 3: In (9), let ( A,B,C ) ASPR, Aq stable, (f) D bounded. Then the Mod- 
ified Adaptive Controller with RMF in (19)-(20) produces e=y and e Q ulti- 


mately bounded into a ball of radius R,=S 



1 Pimx ^ 


Prmn 


M v with exponential 


rate and bounded adaptive gains ( G e ,G D ). 

Proof: Fro m ( 1 9), we have u p = GJe + G D (j) D , s o we can write 
Am = u - m* 

= [G e e y + G D </> D ]-[S* 2 L]</> D 
- G*e y + AG 77 
A G e =G e -G* e 
A G d = G d -(S* 2 L) 

AG = G-G„- [AG e AG d ] . Then 
e„ 


where 


n = 


y 

L'PdA 


with C, 


\ \ 


j C, = A £ + B Au = A c £ + Bw + sT 0 u d 
[e,=CC 

w = AGtj,A c =A+BG* e C . 


From (20), we can see that 

G = AG = -e r/ T h;h = 


K 0 

0 h n 


>0 


(40) 


( 41 ) 
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Since ( A,B,C ) is ASPR, and by the lemma, so is ( A,B,C ), we can we can 
use the following result from [8] where v = T 0 u D is bounded because the distur- 
bance u D —L(j) D is bounded. 

Result: Consider the nonlinear, coupled system of differential equations, 

'g = A c g + B(G{t)-G*)ri+£v 

<e y =Cg (42) 

G(t) = -e y rf h - aG(t) 

_ * 

where G is any constant matrix and h is any positive definite constant matrix, 
each of appropriate dimension. Assume the following: 

i) the triple (A, B ,C) is SPR, 

/ # 'J' * 

ii) there exists Mk> 0 such that IG I G <M K , using the trace norm, 

Hi) there exists M v > 0 such that sup||l/(f)| <M V , 


iv) there exists a > 0 such that a < ■ 


v) h satisfies h 1 < 


2 l aM k 


where p mjn , p max are the minimum and maximum 


eigenvalues of P and q m j„ is the minimum eigenvalue of Q in the system 

\a t c p+pa c =-q 

I PB = C T 

Then the matrix G(t) is bounded and the state tf(t) exponentially approaches the 
ball of radius 

R* = e — — | 7 with s > 0. 


From this result, we have fi, is ultimately bounded into the ball of radius R *, 
which leads to e = y p = y p — y* = CA and e Q ultimately bounded as well. 

Therefore G — G t + A G is bounded, as desired. # 

( 1 + Va^T) 

rt-noormonl-Kr tVio rorline rtf tlio oiror K oil /? = C ^ A jg dCtCf" 


Consequently, the radius of the error ball R»=S- 




mined by the size of e, which is related to the amount of disturbance leakage into 
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the Q modes. It can be seen that, when there is no leakage of the disturbance into 
the Q modes (£’ = 0), the convergence is asymptotic to zero. 

Also, when V = B and Tq =B 0 , it is possible to choose Si =0 and S 2 =— 6 
in (34). Then, even if £ — 1 , the tracking error will asymptotically go to zero. 


Simulation Results with RMF 


In this section we will apply the above theoretical results to a very simple flexible 
structure example to illustrate the behavior with and without the Residual Mode 
Filter. The structure has a rigid body mode and two flexible modes given by: 


1 + 5 


1 


-5 4 +35 3 +05 2 + 35 + l 


-35 3 +25 2 


P(s) = 

5 Z 5 Z +5 + l 5^+5 + 2 5°+25 3 + 45^ 

This example can obviously be extended to have many more flexible modes. But 
we are only trying to illustrate the value of the RMF approach. More seriously 
complex flexible structures are being addressed but will have to await luture pa- 
pers. 

This plant has two non-minimum phase zeros at 0.422+0.9543/ and thus does 


not meet 

w=- 


the 

5 


ASPR condition. However, when the middle mode 


5 +5 + 1 

is removed, the plant becomes: P{s) — 


1 + 5 


-5 + 2 


■3 s +3s+2 

+ 5 3 + 25 2 


which is minimumphase and has a state space realization given by: 


| X p = Ax p 

= C *p 


-B(u p 


>) 


with A — 


0 

1 

0 

0 


0 


2 

0 

0 

1 

0 


0 


3 





, r = 5= 


, c T = 


0 

0 

0 

1 


0 


3 

_0 

0 

-2 

-1. 


_i. 


.1. 


with C2? = 1, so CB is nonsingular. Therefore, ( A,B,C ) is ASPR. 

The reference model to be tracked is 

[<c m =(-\)x m +(\)u m 

\y m = OK 


which is excited by steps generated by U m = (0 )u m . The matching conditions are 
known to be solvable, but their solution is not needed to apply the theory. 
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The RMF generated by P 0 (s) = 


S“+5 + l 


is represented by Aq — 


0 

-1 


1 

-1 



c e = [-3 


0]. And we see that C 0 B 0 — 0. 


The adaptive controller given by (38) - (39) is implemented with /7 U = 10, h m = 1, 
h e = 1 0 , /7 d= 100, and a= 0. The disturbance is a nondimens ional step of size 10. Set- 
ting e=l, we obtain figs. 1 and 2 from a MatLab/Simulink simulation. The output 
tracking error is shown to converge to zero in fig. 1. The adaptive gains also con- 
verge in fig. 2. This illustrates the behavior of the adaptive controller plus the 
second order RMF. Without the RMF, the plant and adaptive controller are imme- 
diately unstable in closed-loop. 


Conclusions 

We have proposed a modified adaptive controller with a residual mode filter. The 
RMF is used to accommodate problematic modes in the system that inhibit the 
adaptive controller, in particular the ASPR condition. This new theory accounts 
for adaptive model tracking and for leakage of the disturbance term into the Q 
modes. A simple three mode example shows that the RMF can restore stability to 
an otherwise unstable adaptively controlled system This is done without seriously 
modifying the adaptive controller design. 
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Adaptive Output Tracking with RMF 
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Fig. 1 . Nondimensional output tracking response with adaptive controller augmented with RMF. 



Fig. 2. Adaptive gains, Ge=error gain, Gd=disturbance gain. 


Appendix 


Proof of Lemma 1 : The Linear Matching Conditions (8) can be rewritten: 

\AS l+ BS 2 = S l L m+ H l 

1 cs x =h 2 
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where S^\S* U S* l2 S^] S 2 = [sT 2i S* 22 S* 22 ] L n 


A m B m 0 

o F m 0 

0 0 F 


and 


//, = [0 0 TO] 

h 2 = \C m 0 0]‘ 

Suppose CB is nonsingular. Use the coordinate transformation W from Lemma 
2 in [11] to put (A, B, C) into normal form: 

I y = A n y + A l2 z 2 + CBu 
|z 2 = A 21 y + A 22 z 2 

i.e., 


WB = 




c 



A n 

A I2 

there exists 

W = 


such that WAW 

l =A = 



L w t 2 p 2 \ 



- A 21 

A 22- 


CB 

0 


= B, and CW 1 =[/,„ 0] = C which implies that 


S 1 L m = WS 1 L m = WAW^WS X + WBS 2 - WH 1 = AS X + BS 2 -H x 

~cb\ 


- AS, 


0 


I s 2 -h x 


and H 2 =CW 'BS, = CS ] =\l 0].S, = S a where 5, ee 

which implies that 




From 



H-, 


H 2 


CB 


H, 

this we have that 

Z 

IaJ 

K= A 

J b \ 

+ 

_ 0 _ 

S 2 - 

Cl 

H b . 


S 1 ={CB)-\H 1 L m +H a ~(A n H 2 +A l2 S b )] 

\s b L m =A 22 SbHA 2l H 2 -H h ) 

Now, if ( A 22 ,L m ) share no eigenvalues, it is well known [5] that we can solve 

c - \ H 2 

the above for a unique o b and conveisely, then S x = _ 

Sb 

S 2 =(CB) \H 2 L m +H a —(A ll H 2 +A l2 S b )] and AS b —S b L m =H b . Since 
( A 22 ,L m ) share no eigenvalues, this is the same as A 22 sharing no eigenvalues 
with A m , F m or F . But the eigenvalues of A 22 from its normal form are known 
to be the transmission zeros of the open -loop system(zf,i?,C) ; see e.g. [13]. Thus, 
we have proved the result. End of Proof. 


Proof of Theorem 2 
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It was already shown that e* and AG are bounded. To prove that 
e * — - — > 0, we must use the following version of Barbalat’s lemma; see [19] 

pp. 210-211: 

Lemma'. If f(t) is a real, differentiable function on (0,co) with lim / ( t) finite 

» 00 

df .. df 

and — uniformly continuous, then lim — = 0. 

dt '-*» dt 

We have already seen that fift) < 0; therefore V(t) — V( 0) = f V(z)dT< 0 

JO 

or 0<V(t)<V(0 ) where F(0) <oo. Hence lim V(t) is finite. Also, V (t) is 

t —> 00 

bounded because 

V(t ) = -(e/Qe*)< ||e»||||£)||||e„|| 


=|e,J|g||||v4 c e. + BAGr/j 
<||e»||||0||(||^c|||| e *|| + ||5|||AG|||7||) 

and e, and AG are bounded by the previous argument via Lyapunov theory. Also 
tj is bounded since u m is bounded, A m is stable, e y —Ce * is bounded, and (j. b D is 

bounded. Thus V(t) = [ V(z)dT is uniformly continuous and Barbalat’s 


Lemma may be applied to yield: 

0 = lim v(t) = — Umle/pe* ). Since (?>0,wehave e* >0, as desired. 

t->oo t —> oo ^ ' t—*x> 

End of Proof. 



